(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 8.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[      8982,        309]
NotebookOptionsPosition[      7895,        267]
NotebookOutlinePosition[      8233,        282]
CellTagsIndexPosition[      8190,        279]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 RowBox[{"N1", "=", 
  RowBox[{"1", "-", "\[Xi]", "-", "\[Eta]"}]}], "\[IndentingNewLine]", 
 RowBox[{"N2", "=", "\[Xi]"}], "\[IndentingNewLine]", 
 RowBox[{"N3", "=", "\[Eta]"}]}], "Input",
 CellChangeTimes->{{3.546090090259152*^9, 3.546090117225739*^9}}],

Cell[BoxData[
 RowBox[{"1", "-", "\[Eta]", "-", "\[Xi]"}]], "Output",
 CellChangeTimes->{3.546090118343816*^9, 3.546090227791521*^9}],

Cell[BoxData["\[Xi]"], "Output",
 CellChangeTimes->{3.546090118343816*^9, 3.546090227796844*^9}],

Cell[BoxData["\[Eta]"], "Output",
 CellChangeTimes->{3.546090118343816*^9, 3.546090227803439*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"\[DoubleStruckX]", "=", 
  RowBox[{
   RowBox[{"N1", " ", "x1"}], "+", 
   RowBox[{"N2", " ", "x2"}], "+", 
   RowBox[{"N3", " ", "x3"}]}]}], "\[IndentingNewLine]", 
 RowBox[{"\[DoubleStruckY]", "=", 
  RowBox[{
   RowBox[{"N1", " ", "y1"}], "+", 
   RowBox[{"N2", " ", "y2"}], "+", 
   RowBox[{"N3", " ", "y3"}]}]}]}], "Input",
 CellChangeTimes->{{3.546090128781515*^9, 3.546090176510607*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"x3", " ", "\[Eta]"}], "+", 
  RowBox[{"x1", " ", 
   RowBox[{"(", 
    RowBox[{"1", "-", "\[Eta]", "-", "\[Xi]"}], ")"}]}], "+", 
  RowBox[{"x2", " ", "\[Xi]"}]}]], "Output",
 CellChangeTimes->{{3.546090152716256*^9, 3.546090178617481*^9}, 
   3.546090229653254*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"y3", " ", "\[Eta]"}], "+", 
  RowBox[{"y1", " ", 
   RowBox[{"(", 
    RowBox[{"1", "-", "\[Eta]", "-", "\[Xi]"}], ")"}]}], "+", 
  RowBox[{"y2", " ", "\[Xi]"}]}]], "Output",
 CellChangeTimes->{{3.546090152716256*^9, 3.546090178617481*^9}, 
   3.5460902296584187`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"sol", "=", 
  RowBox[{
   RowBox[{"Solve", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"\[DoubleStruckX]", "\[Equal]", "x"}], ",", 
       RowBox[{"\[DoubleStruckY]", "\[Equal]", "y"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"\[Xi]", ",", "\[Eta]"}], "}"}]}], "]"}], "[", 
   RowBox[{"[", "1", "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5460901801165*^9, 3.546090245085465*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"\[Xi]", "\[Rule]", 
    RowBox[{"-", 
     FractionBox[
      RowBox[{
       RowBox[{"x1", " ", "y"}], "-", 
       RowBox[{"x3", " ", "y"}], "-", 
       RowBox[{"x", " ", "y1"}], "+", 
       RowBox[{"x3", " ", "y1"}], "+", 
       RowBox[{"x", " ", "y3"}], "-", 
       RowBox[{"x1", " ", "y3"}]}], 
      RowBox[{
       RowBox[{"x2", " ", "y1"}], "-", 
       RowBox[{"x3", " ", "y1"}], "-", 
       RowBox[{"x1", " ", "y2"}], "+", 
       RowBox[{"x3", " ", "y2"}], "+", 
       RowBox[{"x1", " ", "y3"}], "-", 
       RowBox[{"x2", " ", "y3"}]}]]}]}], ",", 
   RowBox[{"\[Eta]", "\[Rule]", 
    RowBox[{"-", 
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"-", "x1"}], " ", "y"}], "+", 
       RowBox[{"x2", " ", "y"}], "+", 
       RowBox[{"x", " ", "y1"}], "-", 
       RowBox[{"x2", " ", "y1"}], "-", 
       RowBox[{"x", " ", "y2"}], "+", 
       RowBox[{"x1", " ", "y2"}]}], 
      RowBox[{
       RowBox[{"x2", " ", "y1"}], "-", 
       RowBox[{"x3", " ", "y1"}], "-", 
       RowBox[{"x1", " ", "y2"}], "+", 
       RowBox[{"x3", " ", "y2"}], "+", 
       RowBox[{"x1", " ", "y3"}], "-", 
       RowBox[{"x2", " ", "y3"}]}]]}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.546090213841282*^9, 3.546090245563058*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"\[DoubleStruckCapitalN]1", "=", 
  RowBox[{
   RowBox[{"N1", "/.", "sol"}], "//", "Simplify"}]}], "\[IndentingNewLine]", 
 RowBox[{"\[DoubleStruckCapitalN]2", "=", 
  RowBox[{
   RowBox[{"N2", "/.", "sol"}], "//", "Simplify"}]}], "\[IndentingNewLine]", 
 RowBox[{"\[DoubleStruckCapitalN]3", "=", 
  RowBox[{
   RowBox[{"N3", "/.", "sol"}], "//", "Simplify"}]}]}], "Input",
 CellChangeTimes->{{3.5460902478761168`*^9, 3.546090292178389*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"x2", " ", "y"}], "-", 
   RowBox[{"x3", " ", "y"}], "-", 
   RowBox[{"x", " ", "y2"}], "+", 
   RowBox[{"x3", " ", "y2"}], "+", 
   RowBox[{"x", " ", "y3"}], "-", 
   RowBox[{"x2", " ", "y3"}]}], 
  RowBox[{
   RowBox[{"x2", " ", "y1"}], "-", 
   RowBox[{"x3", " ", "y1"}], "-", 
   RowBox[{"x1", " ", "y2"}], "+", 
   RowBox[{"x3", " ", "y2"}], "+", 
   RowBox[{"x1", " ", "y3"}], "-", 
   RowBox[{"x2", " ", "y3"}]}]]], "Output",
 CellChangeTimes->{{3.546090253521658*^9, 3.546090292784772*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"x3", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "y"}], "+", "y1"}], ")"}]}], "+", 
   RowBox[{"x1", " ", 
    RowBox[{"(", 
     RowBox[{"y", "-", "y3"}], ")"}]}], "+", 
   RowBox[{"x", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "y1"}], "+", "y3"}], ")"}]}]}], 
  RowBox[{
   RowBox[{"x3", " ", 
    RowBox[{"(", 
     RowBox[{"y1", "-", "y2"}], ")"}]}], "+", 
   RowBox[{"x1", " ", 
    RowBox[{"(", 
     RowBox[{"y2", "-", "y3"}], ")"}]}], "+", 
   RowBox[{"x2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "y1"}], "+", "y3"}], ")"}]}]}]]], "Output",
 CellChangeTimes->{{3.546090253521658*^9, 3.5460902927897863`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"x1", " ", "y"}], "-", 
   RowBox[{"x2", " ", "y"}], "-", 
   RowBox[{"x", " ", "y1"}], "+", 
   RowBox[{"x2", " ", "y1"}], "+", 
   RowBox[{"x", " ", "y2"}], "-", 
   RowBox[{"x1", " ", "y2"}]}], 
  RowBox[{
   RowBox[{"x2", " ", "y1"}], "-", 
   RowBox[{"x3", " ", "y1"}], "-", 
   RowBox[{"x1", " ", "y2"}], "+", 
   RowBox[{"x3", " ", "y2"}], "+", 
   RowBox[{"x1", " ", "y3"}], "-", 
   RowBox[{"x2", " ", "y3"}]}]]], "Output",
 CellChangeTimes->{{3.546090253521658*^9, 3.546090292795383*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"D", "[", 
   RowBox[{"\[DoubleStruckCapitalN]1", ",", "x"}], "]"}], "//", 
  "Simplify"}]], "Input",
 CellChangeTimes->{{3.546090295279338*^9, 3.546090314109045*^9}, {
  3.546090408171804*^9, 3.546090408381731*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"-", "y2"}], "+", "y3"}], 
  RowBox[{
   RowBox[{"x3", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "y1"}], "+", "y2"}], ")"}]}], "+", 
   RowBox[{"x2", " ", 
    RowBox[{"(", 
     RowBox[{"y1", "-", "y3"}], ")"}]}], "+", 
   RowBox[{"x1", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "y2"}], "+", "y3"}], ")"}]}]}]]], "Output",
 CellChangeTimes->{{3.546090301839141*^9, 3.546090314721261*^9}, {
  3.546090404759035*^9, 3.5460904089251757`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"D", "[", 
   RowBox[{"\[DoubleStruckCapitalN]1", ",", "y"}], "]"}], "//", 
  "Simplify"}]], "Input",
 CellChangeTimes->{{3.54609041321425*^9, 3.546090420106884*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{"x2", "-", "x3"}], 
  RowBox[{
   RowBox[{"x2", " ", "y1"}], "-", 
   RowBox[{"x3", " ", "y1"}], "-", 
   RowBox[{"x1", " ", "y2"}], "+", 
   RowBox[{"x3", " ", "y2"}], "+", 
   RowBox[{"x1", " ", "y3"}], "-", 
   RowBox[{"x2", " ", "y3"}]}]]], "Output",
 CellChangeTimes->{{3.546090410626383*^9, 3.546090420594162*^9}}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.546090219753129*^9, 3.546090224122991*^9}, 
   3.54609031824605*^9}]
},
WindowSize->{740, 623},
WindowMargins->{{142, Automatic}, {Automatic, 24}},
FrontEndVersion->"8.0 for Linux x86 (64-bit) (February 23, 2011)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 271, 5, 69, "Input"],
Cell[853, 29, 133, 2, 30, "Output"],
Cell[989, 33, 96, 1, 30, "Output"],
Cell[1088, 36, 97, 1, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[1222, 42, 419, 11, 50, "Input"],
Cell[1644, 55, 301, 8, 30, "Output"],
Cell[1948, 65, 303, 8, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2288, 78, 440, 12, 30, "Input"],
Cell[2731, 92, 1300, 38, 87, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4068, 135, 465, 10, 69, "Input"],
Cell[4536, 147, 548, 16, 49, "Output"],
Cell[5087, 165, 713, 25, 50, "Output"],
Cell[5803, 192, 548, 16, 49, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6388, 213, 250, 6, 30, "Input"],
Cell[6641, 221, 520, 17, 50, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7198, 243, 200, 5, 30, "Input"],
Cell[7401, 250, 358, 10, 49, "Output"]
}, Open  ]],
Cell[7774, 263, 117, 2, 30, "Input"]
}
]
*)

(* End of internal cache information *)
